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Abstract 

The uncertainty in the contribution to the Lamb shift in muonic hydrogen, AE subt arising from 
proton polarizability effects in the two-photon exchange diagram at large virtual photon momenta 
is shown large enough to account for the proton radius puzzle. This is because AE subt is deter- 
mined by an integrand that falls very slowly with very large virtual photon momenta. We evalu- 
ate the necessary integral using a set of chosen form factors and also a dimensional regularization 
procedure which makes explicit the need for a low energy constant. The consequences of our 
two-photon exchange interaction for low-energy elastic lepton-proton scattering are evaluated 
and could be observable in a planned low energy lepton-proton scattering experiment planned to 
run at PSI. 



I. INTRODUCTION 

The proton radius puzzle is one of the most perplexing physics issues of recent times. 
The extremely precise extraction of the proton radius [T] from the measured energy dif- 
ference between the 2P0^ 2 and 2Sj^ states of muonic hydrogen disagrees with that ex- 
tracted from electronic hydrogen. The extracted value of the proton radius is smaller than 
the COD ATA Q value (based mainly on electronic H) by about 4% or 5.0 standard devi- 
ations. This implies [1J that either the Rydberg constant has to be shifted by 4.9 standard 
deviations or that present QED calculations for hydrogen are insufficient. The Rydberg 
constant is extremely well measured and the QED calculations seem to be very extensive 
and highly accurate, so the muonic H finding is a significant puzzle for the entire physics 
community. 

Pohl et al. show that the energy difference between the 2P^ 2 and 2S[^ states, AE is 
given by 

AE = 209.9779(49) - 5.2262r 2 + 0.0347rJ meV, (1) 

where r-p is given in units of fm. Using this equation and the experimentally measured 
value AE = 206.2949 meV, one can see that the difference between the Pohl and COD ATA 
values of the proton radius would be removed by an increase of the first term on the rhs 
of Eq. (1) by 0.31 meV=3.1 x 10" 10 MeV. 

This proton radius puzzle has been attacked from many different directions |3|-||2TH 
The present communication is intended to investigate the hypothesis that the proton po- 
larizability contributions entering in the two-photon exchange term, see Fig. [TJ can ac- 
count for the 0.31 meV. This idea is worthy of consideration because the computed effect 
is proportional to the lepton mass to the fourth power, and so is capable of being relevant 




FIG. 1: The box diagram for the 0{a?m i ) corrections. The graph in which the photons cross is 
also included. 
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for muonic atoms, but irrelevant for electronic atoms. 

II. AE subt AND ITS EVALUATION 

The basic idea is that the two-photon exchange term depends on the forward virtual 
Compton scattering amplitude T^ v (v,q 2 ) where q 2 is the square of the four momentum, 
ql 1 of the virtual photon and v is its time component. One uses symmetries to decom- 
pose Tt lv (v,q 2 ), into a linear combination of two terms, 7^2 (v,^ 2 ). The imaginary parts 
of Ti ^(v,^ 2 ) are related to structure functions Fx 2 measured in electron- or muon-proton 
scattering, so that T\ % can be expressed in terms of F1.2 through dispersion relations. 
However, V\{y, Q 2 ) falls off too slowly for large values of v for the dispersion relation to 
converge. Hence, one makes a subtraction at v = 0, requiring that an additional function 
of Q 2 (the subtraction function) be introduced. One accounts for the nucleon Born terms, 
and the remainder of the unknown subtraction function is written as T\ (0, Q 2 ) |fl7f. This 
term is handled by making a power series expansion around Q 2 = 0, and then using 
effective field theory to determine the coefficients of the series. The problem with using 
this expansion is that this contribution to the energy is determined by an integral over all 
values of Q 2 . We proceed by elaborating the consequences of the behavior of Ti(0, Q 2 ) 
for large values of Q 2 . This is followed by the development of an alternate effective field 
theory approach to the muon-proton scattering amplitude. In either case, one can ac- 
count for the needed Lamb shift, while also providing consequences for the two-photon 
exchange contribution to the scattering amplitude that can be tested in an upcoming ex- 
periment 1122 II . 

The contribution to the Lamb shift that is caused by T\ (0, Q 2 ) is denoted as AE subt and 
is given by |HzHl3|23r|23 

AE^ = ^m[%K^)M^Q 2 ) (2) 

3 3 

where (p 2 (0) = for the 2S state with m, (m r ) as the lepton (reduced) mass, and 

h{t) = (1 - It) ((1 + -) 1/2 - l) + 1. (3) 

The function h(t) is monotonically falling, approaching 1/ \ft for small values of t, and 
falling as 3/ (4f) large values of t. The subtraction function Ti(0, Q 2 ) is not available 

3 



from experimental measurements, except at the real photon point Q 2 = 0. It comes from 
the excitation of the proton, and can be described, at small values of Q 2 , in terms of the 
electric (oce) and magnetic (/3m) polarizabilities. For small values of Q and v = one 
sees [23] lim y 2 q2^ Ti(0, Q 2 ) = ^-/$m- Using this simple linear Q 2 -dependence inEq. (2) 
shows that the integral over Ti(0, Q 2 ) converges at the lower limit, but diverges loga- 
rithmically at the upper limit. Thus obtaining a non-infinite result depends on including 
an arbitrary form factor that cuts off the integrand for large values of Q 2 or some other 
renormalization procedure. 

We note that limg2_ > . 00 Ti(0, Q 2 ) can be obtained from the operator production expan- 
sion |26]|27]|. Using Eq. (2.18) of Ref. [26], neglecting the term proportional to light quark 
masses, and accounting for different conventions yields Ti(0, Q 2 ) ~ 2.1 fm~V Q 2 . This 
1/Q 2 behavior removes the putative logarithmic divergence of Ti(0, Q 2 ), but this func- 
tion is far from determined. 

We follow the previous literature by including a form factor defined as Fi oop . Then 

T 1 (0,Q 2 ) = ^Q 2 F loop (Q 2 ). (4) 
Using Eqs. ( 2|3|4 > one finds the energy shift to be 



AE suU = ft¥(0) gM ; dQf 
m ct 



:i-2Q 2 /(4m 2 )) h/l + ^-l| +1 



fioo P (Q 2 )- 

(5) 



The issue here is the arbitrary nature of the function Fi 00 p(Q 2 )- Pachucki ||2"4H used the 
dipole form, ~ 1/ Q 4 , often used to characterize the proton electromagnetic form factors. 
But the subtraction function should not be computed from the proton form factors, be- 
cause virtual component scattering includes a term in which the photon is absorbed and 
emitted from the same quark [28]. Carlson and Vanderhaeghen [fT7| evaluated a loop di- 
agram using a specific model and found a form factor ~ 1 / Q 2 log Q 2 , leading to a larger 
contribution to the subtraction term than previous authors. Birse & McGovern [20] use 
terms up to fourth-order in chiral perturbation theory to find 

Tr(o,Q 2 )-^Q 2 fi-^ + 0(Q 4 )) -^Q 2 1 x „ (6) 
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with Mp = 460 ± 50 MeV. They also use the most recent evaluation of /3 m, based on a fit 
to real Compton scattering 1129 II that finds 

f> M = (3.1 ± 0.5) x 10~ 4 fm 3 , (7) 

where only statistical and Baldin Sum Rule errors are included. Their result is a negligible 
AE subt = 4.1/^ eV Il20l . The form Eq. ^ achieves the correct 1/Q 2 asymptotic behavior 
of Ti(0, Q 2 ) but the coefficient Pm/k is n °t the same as obtained from the operator prod- 
uct expansion. The coefficient of Eq. (|6j) is about twice the asymptotic limit obtained by 
Collins EU. 

Previous authors fT7ll20| noted the sensitivity of the integrand of Eq. ^ to large values 
of Q 2 . Our aim here is to more fully explore the uncertainty in the subtraction term 
that arises from the logarithmic divergence. We shall use a form of fi op(Q 2 ) that is 
consistent with the constraint on the Q 4 term found Birse & McGovern II20II . This is done 
by postulating a term that begins at order Q 6 in Eq. (jij), such as 

(2 \ n 
%A (TTW'"- 2 ' N -" + 3 ' (8) 



where Mo, a are parameters to be determined. With Eq. ^ the low Q 2 behavior of 
Ti (0, Q 2 ) is of order Q 6 or greater and it falls as 1 / Q 4 or greater for large values of Q 2 . So 
far as we know, there are no constraints on the coefficient of the Q 6 term and the 1/Q 4 
term. However, we shall determine the subtraction term's contribution to the Lamb shift 
as a general function of n, N. We note that f$M is anomalously small due to a cancellation 
between pion cloud and intermediate A terms p0| , so that one can use a value ten times 
larger than appears in Eq. (j7j) to set the overall scale of the subtraction term. Thus we 
replace the term of Eq. Q by a general form of the same dimensions /3: fi^ — > ft- 

The use of Eq. (|8|) in Eq. ^ allows one to state the expression for the energy shift in 
closed form as a general function of n, N. We find 

iP-^jf')",^), (9) 



h,,N(m 2 a) , if * ^ [(1 - ^) (d + ^ -l) +1 



(10) 



The integral over x can be obtained in a closed form in terms of hypergeometric func- 
tions. However, a much more understandable expression can be obtained by replacing 



the bracketed expression in Eq. (10 1 by its large argument limit (3m 2 a/x). This approxi- 
mation is valid over the entire range of the integrand because of the presence of the factor 
x n with n > 2. Then one obtains 

h, N {m 2 a) « 3m 2 r(M r " ( ^ r(n) = 3m 2 B(N,n), (11) 

so that 



AE subt « 3oc 2 m(p 2 (0) —A n B(N,n), A = (12) 



Numerical evaluations show that the approximation is accurate to better than a quarter 
of a percent. The expression Eq. ( 12 1 makes clear the m 4 dependence of the contribution 



to the Lamb shift. 

The numerical value of the term AE subt depends on (n,N),fi and the combination 

Mia = A" 1 : 

AE = 3.91meV fm 3 J 6A"B(N,n). (13) 

If we take N = 5,n = 2 so that B(5,2) = 1/12, and 6 = 10" 3 fm" 3 , a value of A = 30.9 
reproduces E = 0.31 meV. If we take M = 0.5 GeV (as in flU) , then a" 1 = 15.4 GeV 2 , 
and that the contribution to the integral comes from the region of very high values of Q 2 . 
Other values of n, N and A could be used to get the identical contribution to the Lamb 
shift. 

Chiral perturbation theory could be used to determine the terms of order Q 6 and 
higher in rf M (0, Q 2 ), but this procedure is always limited to a finite number of terms. 
Indeed one could use values of n greater than 2, and still reproduce the needed contribu- 
tion to the Lamb shift. 

The above discussion shows that the current procedure used to estimate the size of the 
subtraction term is rather arbitrary. This arises because the chiral EFT is being applied 
to the virtual-photon nucleon scattering amplitude. Another technique would be to de- 
velop an effective field theory to determine the short-distance lepton-nucleon amplitude 
implied by the subtraction term. 
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III. EFFECTIVE FIELD THEORY FOR THE ]if INTERACTION 



The previous considerations show that the value of AE suht depends heavily on as- 
sumptions the behavior of 7i(0, Q 2 ), (F| rm / oop ) for large values Q 2 . This is true even 
though the leading 1/Q 2 term is known. The underlying cause of this uncertainty is the 
would-be logarithmic divergence in the integral of Eq. (|5j) for the case Fi oop = 1. This is 
a symptom that some other technique could be used |3~T|. Another way to proceed is to 
use an effective field theory (EFT) for the lepton-proton interaction [32 J. In EFT, logarith- 
mic divergences identified through dimensional regularization are renormalized away 
by including a lepton-proton contact interaction in the Lagrangian. 

We may handle the divergence using standard dimensional regularization (DR) tech- 
niques by evaluating the scattering amplitude of Fig. 1. The term of interest is obtained 
by including only T\{Q, Q 2 ) of Eq. |3j) with Fi oop = 1. We evaluate the loop integral in 
d = 4 — e dimensions and obtain the result: 

M£ R (loop) = \i a 2 m^ [- + log ^ + % - 7e + \o%An\u fUi U f U u (14) 
2 ol L e m z 6 J 1 1 

where lower case spinors represent leptons of mass m, and upper case proton of mass M, 

q is momentum transferred to the proton, and 7e is Euler's constant, 0.577216- • • . 

The result Eq. ( [14] > corresponds to an infinite contribution to the Lamb shift in the limit 

that e goes to zero. In EFT one removes the divergent piece by adding a lepton-proton 

contact interaction to the Lagrangian that removes the divergence, replacing it by an 

unknown finite part. The finite part is obtained by fitting to a relevant piece of data. Here 

the only relevant data is the 0.31 meV needed to account for the proton radius puzzle. 

The low energy term contributes 

A4f R (LET) = iC(}i), (15) 



where C(ji) is chosen such that the sum of the terms of Eq. (14) and Eq. (15), = Ai^' i s 
finite and independent of the value of ]i. Thus we write the resulting scattering amplitude 
as 

M% R = i a 2 m^-{X + 5/4) u f UiU f Ui, (16) 



where A is determined by fitting to the Lamb shift. Eq. ( 16 1 corresponds to using the MS 



scheme because the term log(47r) — 7£ is absorbed into A. 



The corresponding contribution to the Lamb shift is given by 

AE DR = a 2 m^ 2 (0)(A + 5/4). (17) 

Setting AE DR to 0.31 meV in the above equation requires that A = 769, which seems like 
a large number. However, as noted above, ($m is extraordinarily small II30II . The natural 
units of polarizability are ^ ~ 4n/ A^, 1133 II where A^ = 4nf n , (f n is the pion decay 
constant). Then Eq. ( [l6| > becomes 



M% R = i 3.95 a 2 m-^u f UiU f Ui. (18) 



The coefficient 3.95 is of natural size. Thus standard EFT techniques result in an effective 
lepton-proton interaction of natural size that is proportional to the lepton mass. The form 



of Eq. (18 1 is not unique. There are other possible operators that reduce to that form in 
the low-energy, low-momentum regime of relevance here. 



The present results, Eq. ( 12 1 and Eq. ( |T7| > represent an assumption that there is a lepton- 
proton interaction of standard-model origin, caused by the high-momentum behavior of 
the virtual scattering amplitude, that is sufficiently large to account for the proton radius 
puzzle. Fortunately, our hypothesis can be tested in an upcoming low-energy ]i ± p,e ± p 
scattering experiment | |22 || planned to occur at PSI. 



IV. LEPTON PROTON SCATTERING AT LOW ENERGIES 

Our aim is to determine the consequences of the particular two-photon exchange term 
for lepton-proton scattering at low energies. Our previous attempt [21J implied very large 
corrections to quasi-elastic electron-nucleus scattering that are severely in disagreement 
with experiment p4|. It is necessary to check that a similar large unwanted contribution 
does not appear here. Thus we provide a prediction for the PSI experiment. It is well- 
known that two-photon exchange effects in electron-proton scattering are small at low 
energies. Our contact interaction is proportional to the lepton mass, so it could provide 
a measurable effect for muon-proton scattering but be ignorable for electron-proton scat- 
tering. We shall investigate the two consequences of using form factors (FF) and effective 
field theory (DR). 
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The invariant amplitude is given as Ai fi = M}} + M-f} where the superscripts de- 



note the number of photons exchanged. The first term is given by 

„2 



M 



(i) 



^ 2 i . n u f j v UiU f T^Ui, 



(19) 



ft " " r > 2 + i0' 

where Uj f represent leptons of mass m, Ujf represent the proton of mass M, and q is mo- 
mentum transferred to the proton. The minus sign holds for negatively charged muons, 
and the plus sign for positively charged muons. 



Tt l = 7 ^F 1 (q 2 ) + i—q v F 2 (q 2 ) f 



(20) 



with Pf = Pj + q. The second-order term that arising from the use of Eq. (|4j) and Eq. |8j) , 
which uses form factors (FF) is given by 



2n) 



where 



L 



1 Tj(0,-/c 2 ) 
k 2 + z'O (k - q) 1 + z'O 

2mk 2 + Api • k fi 



UfLiijUfUi, 



(21) 



-4(p r A:) 2 + fc 4 + i0" 

The second-order amplitude arising from the use of EFT is given above in Eq. (16). 



(22) 



The cross section depends upon the average over initial and sum over final fermion 
spins, as denoted by an over line. We first obtain 



M 



(i) 



fi 



Standard text-book expressions 
use the Mott form, obtained by ignoring the lepton mass. This is not a good approxima- 
tion for the muons of the experiment [22J which have momenta ranging between 100 and 



200 MeV/ c, and our terms of Eq. (21 ) and Eq. (16) would vanish in that approximation. 
We find 



M 



(i) 
fi 



[l6M\ee> + q 2 /4)(F z - ^F 2 ) + 40^/2 + qV)] ( ~ 2 



2„„2\ 



q 2 + z'O 



(23) 



where 8 is the laboratory scattering angle, and e(e') is the incident (final) lepton labora- 
tory total energy. The present interference term of interest A = 2Real [{■M-r i )*(M.\ i )] is 
obtained using standard trace algebra. We find 



A(FF) = 8MG E (q : 



=Fz<r 



d^k 



i r^-k 2 ) 



q 2 + zO (2tt) 4 J k 2 + z'O (k - q) 2 + z'O (-4(p f • k) 2 + fc 4 + z'0) 



x [2m z k z (4eM + q z ) + 4pi-kk- (P t + P f )q z /2 + 2p { ■ kk • (p f + p f )(4eM + tf)). (24) 

9 




FIG. 2: The ratio R obtained using form factor (FF) regularization or dimensional regularization 
(DR). The solid curves show the results for muon laboratory momentum of 200 MeV/c and the 
dashed curves show the results for 100 MeV/ c. 



We have seen that the integrand is dominated by large values of k, therefore we neglect q 
in the integrand. This allows considerable simplification so that we find 



e 2 

A(FF) = T96G E (q 2 )M 2 ^ — -eccBm 2 X n B{N ,n) (25) 



where the terms A, B(N, n) appear in Eq. (12 1. A negligible term proportional to the 
square of the incident lepton momentum has been dropped. The term A adds to the 
square of the lowest order term for fi~p interactions, as expected from an attractive inter- 
action that increases the Lamb shift. The computed value of A(FF) does not depend on 



n, N, A for those values that reproduce the needed Lamb shift via Eq. ( 12 1. 

For EFT the contribution to the cross section via interference can be worked out using 



10 



Eq. flgk to be 



A DR = T 8[4eM + q z ]oc(A + j)m z p M G E {q z )M- 



4 y "~ r™-*v-, '-'Y + iO" 

We are now prepared to display the effects of our two-photon exchange term on p 
scattering at low energies. The size of the effect is represented by the ratio R, with 



Mf 



(26) 
(27) 



The ratio R > for p — p scattering. The numerator of Eq. (j27j) is obtained from either 
Eq. (25 1 (FF) or Eq. (26 1 (DR). The ratio R is proportional to the square of the lepton mass, 
which is negligible for — p scattering. We consider two muon momenta 100 and 200 
MeV/ c. The results are shown in Fig. |2j The angular dependence is dominated by the 
Q 2 = — q 2 term inherent in Eq. (27). The two sets of curves are very similar because the 
size of the effect is constrained by the required energy shift of 0.31 meV. The size of the 
effect should be detectable within the expected sub-1 % accuracy of the PSI experiment. 
We emphasize that our calculation is valid only at low muon laboratory energies. 



V. SUMMARY AND DISCUSSION 



The findings of this paper can be summarized with a few statements: 

• The integrand (see Eq. (j2|) that determines the value of AE subt values slowly with 
large values of Q 2 , causing the uncertainty in the evaluation to be large enough to 
account for the proton radius puzzle. 

• The integrand can be evaluated using one of an infinite set of possible form factors 
or a dimensional regularization procedure. 

• Either method can be used to account for the proton radius puzzle and predict an 
observable effect of a few percent for low energy p — p scattering. 

The literature |T | -||2T | poses several explanations for the proton radius puzzle: The 
electronic-hydrogen experiments might not be as accurate as previously reported, p. — e 
universality might be violated, and that a strong interaction effect entering in a loop dia- 
gram is important for muonic hydrogen, but not for electronic hydrogen. It is beyond the 
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scope of the present paper to argue for the unique correctness of any one of these ideas. 
The strong-interaction effect discussed here is large enough to be testable experimentally, 
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